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Abstract 

The discrepancy on ''Li and ^Li abundances between the observational data and the standard Big 
Bang Nucleosynthesis theory prediction has been a nagging problem in astrophysics and cosmology, 
given the highly attractive and succesful Big Bang paradigm. One possible solution of this lithium 
problem is through hadronic decays of a massive metastable particle which alter the primordial 
element abundances. We explore this possibility using gravitino dark matter framework in which 
the next lightest supersymmetric particle (NLSP) is typically long-lived. We found that stop NLSP 
can provide an attractive solution to the lithium problem. 
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I. INTRODUCTION 



The identity of dark matter is one of the most important questions we currently have in 
astrophysics and cosmology. There are many theoretical candidates for dark matter particles 
proposed by many beyond the standard model hypotheses in particle physics, and gravitino 
is one of them within the supergravity framework Due to the smallness of the 

coupling between gravitino and matter in supergravity models, ~ 1/Mpi, the next lightest 
supersymmetric particle (NLSP) is typically longlived. This has direct implication on the 
light element abundances if the NLSP decayed around or after the time of the Big Bang 
Nucleosynthesis (BBN). Thus BBN provides a stringent constraint on gravitino dark matter 



scenarios 
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HQ . However, on the other hand it seems that the 
standard BBN theory cannot explain all of the primordial light element abundances from 
observations, namely the lithium problem described below. This leads to an interesting 
possibility that some non-standard processes were involved in the BBN. In this paper we 
study feasible solution of the lithium problem by gravitino dark matter with stop NLSP 
scenario. 

There are many possibilities for the NLSP with the gravitino LSP, and each has its own 



henomenological signatures. The ones that have been studied are neutralino NLSP 
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13|, stau NLSP |4, Id, ISl, llO|, IjJ, ll3( , stop NLSP [u] and sneutrino NLSP 



4 

131, 
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We do not discuss the details of those phenomenologies here, but concentrate only on the 
BBN effects by looking at the NLSP decays, lifetime and density before the decay. We 
found that stop NLSP is the most interesting scenario in the sense that it provides possible 
solution to the lithium problem. Therefore we focus on stop particle in this paper. 

It has been recognized that the theoretical prediction of ^Li abundance in the standard 
BBN does not agree with the observational results when we adopt the baryon-to-photon ratio 
of WMAP 5-year, r/ = (6.225 ± 0.170) x 10"^° [3]. Using this ratio, the theoretical value 
for ^Li is much larger than the observational one even if we adopt a relatively high value of 



the observational abundance, logiQ(''Li/H) 



obs 



-9.36 ±0.06 [16] ^ Quite recently Refs [1 



reported that the discrepancy gets worse if we adopt an updated theoretical calculation 
for the reaction rate of ^He(^He,7) ''Be. As for ^Li abundance, on the other hand, recent 



^ On the other hand in Ref. [l7| a relatively low value of ^Li abundance was reported, log2o(^Li/H)obs 
—9.90 ± 0.06. Obviously it would be more difficult and problematic to fit this value. 
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observation shows that the theoretical value is much smaller than that of the observation, 



^Li/^Li) 



obs 



0.046 ± 0.022 [19|. We collectively call these discrepancies of Li and Li as 



the "lithium problem" . 

If there is a metastable massive particle which decays, producing energetic standard 
model particles, around or after the BBN era, the light element abundances might be altered 
through electromagnetic and hadronic shower effects. In the hadronic-decay scenario, we 
might be able to solve the lithium problem because the hadron-emission could possibly 
reduce ^Li [2^ and produce ^Li simultaneously 
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22 1 as will be discussed later. To solve 



the lithium problem, the NLSP abundance times its net visible energy which fragments into 
hadrons -Ehad should be in the range of i^had ^nlsp/s — 10^^^ GeV - 10^^^ GeV with their 
lifetime 10^ - 10^ sec 20|, |2l|, |2^ . It is attractive that the stop abundance is naturally tuned 



to solve the lithium problem when we consider a further (second) annihilation which must 
occur just after the QCD phase transition as was pointed out by [23|. ^ 

The outline of this paper is as follows. In section 2 we discuss supersymmetric models 
with stop NLSP, assuming gravitino as the LSP. In section 3 we calculate the stop abundance 
before it eventually decays to gravitino. In section 4 we briefly review the Big Bang Nucle- 
osynthesis theory and then presented our analysis with stop NLSP. Finally, we conclude in 
section 5. 



II. MODELS WITH STOP NLSP 

The lightest stop particle, ti, can be light by (little) seesaw mechanism between left and 
right stop, tL and tR respectively. Recall the stop mass matrix 



ML ML 
ML MIr 



where 



ML 



mIr 

ML 



Mi+ml+^-{Aml. 



m|) cos 2(3 



M~ + -m| cos 2f3 sin^ Ow 

^ 3 

—mt{At + /i cot P) 



(1) 

(2) 

(3) 
(4) 



^ For another simultaneous solutions to solve the lithium problem in astrophysics and cosmology, see also 
and 



Refs. 
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28|, respectively. 
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The off diagonal terms are multiplied by the top-quark mass rrit which is quite large. If At is 
also large, it is possible to have ti lighter than all the other MSSM sparticles, and therefore 
becomes the NLSP (with gravitino as the LSP). 

As we mentioned, in this gravitino LSP scenario, the stop NLSP has a long lifetime. In 
colliders, any metastable stop produced would quickly hadronize either into some sbaryons 
or mesinos. The super-hadrons then decay to the lightest sbaryon tud, which is charged, or 



the lightest mesino, tu, which is neutral 29|] ^. They would then escape the calorimeter, and 
the charged ones could be detected by the muon detector. This hadronization complicates 
the determination of the detection rates of stop. Nevertheless, there has been analysis by 
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CDF at Tevatron which set the lower bound of (meta)stabe stop mass at ~ 250 GeV |3Cl|. 

Realization of stop NLSP scenario in some specific MSSM models was studied in 11 | 
However, it was found that in the Constrained MSSM (CMSSM) it is not possible to have 
stop NLSP due to the Higgs mass and the stop mass lower bound constraints. In the 
Non Universal Higgs Masses (NUHM) model 33|, |3J] we can still have a narrow allowed 



region with stop NLSP. Nonetheless, stop NLSP is still possible if we forgo the universality 



assumption for sfermion and gaugino masses. For example, in [35| it was shown that stop 
would be relatively light if the gluino mass at the GUT scale is lower than the bino and wino 
mass. In that paper the author still assumed neutralino be lighter than the stop. However, 
it would be easy to check that we can combine this non-universality with moderate value of 
trilinear coupling At to get stop be the lightest. Our approach here is not to look at specific 
model of supersymmetry. We will just assume that stop is the NLSP and treat the stop 
mass as a free parameter. 

We then calculate the stop lifetime as a function of stop and gravitino masses. When the 
mass gap between stop and gravitino is larger than rrit the dominant decay channel is the 
2-body decay t G + t. The stop lifetime can then be determined as 



Tj- ~ 487rMpm ~m~ [m~ — m ~ — + 4 sin 6'^cos O^jmtm^^ ^ 

X ^(m~ + — TUty — Am~my\ ^ [{m~ + — m'^)'^ — Am-rrit] ^^"^ , (5) 



For simplicity we will ignore the label '1' to denote the lightest state for the rest of this paper. 
^ There are also studies for future detectabilities of a relatively long-lived stop in the Large Hadron Collider 



(LHC) [3]| and the International Linear Collider (ILC) [32|. However, these studies assume neutralino 
LSP and the analyses were based on the stop decay to charm-quark plus neutralino, and therefore not 
applicable to our case here. 
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FIG. 1: Stop lifetime in the plane of stop mass mj-and gravitino mass m-^i2 = m^. The shadowed 
region on the left with respect to the dotted line is excluded by the experimental lower bound on 
the stop mass rrij > 250 GeV. The left side of the thick line is not allowed by kinematic for 2-body 
decay, t — > t + G. 



where Mp is the reduced Planck mass (Mp = Mpi/y/Sn ~ 2.4 x 10^^ GeV) and % as the 



stop mixing angle. As shown in ll|], the dependence on 6j is small, hence we can fix 6j for 
our analysis (=1.3 radian). Our numerical results are shown in Fig. [H as contours in the mj 
vs mg plane. We do not need to calculate the 3-body decay rate since it is only important 
for lifetime longer than ~ 10^ s. 



III. STOP ABUNDANCE 



We can calculate the relic abundance of the thermally-produced stop according to the 
standard method for calculating the freeze-out value HI]. Assuming that coannihilation 
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effect is negligible and that there is no resonance, we can simplify the calculation by con 



sidering only annihilation channels through gluon exchange. See [36|] for detail discussion 
on this approximation. In this case, the relic density calculation depends only on the stop 
mass. We get a relic density of Q-jrh"^ = 10^^ — 10^^ for = 10^ — 10'^ GeV, corresponding 
to mjY^ = 10~^^ — 10^^^ GeV. However, in addition to the standard annihilation process, 
there is a further annihilation that reduces the number density of stop before it eventually 
decays j^ ]. 

After the QCD phase transition which occurs at a temperature T = Tqcd ~ 150 MeV, 
all of the strong interacting particles are confined into hadrons. The strength of the strong 
interaction is then determined by the effective theory of hadron physics. According to the 
heavy-quark effective theory S^], the length scale for the interactions of heavy-quark hadrons 
is -Rhad ~ 1 - 10 GeV~^ [23|, and we can parameterize the length scale as i?had = fa/^ir 
with is the pion mass ( ~ 140 MeV) and is a numerical parameter ranging from about 
0.1 to 1. 

The relic stop and antistop would also form hadrons, called super-hadrons or shadrons, 
which then undergo further annihilation via bound states. The annihilation rate among 
these shadrons can be written as 

Tann = ipv) 71^ , (6) 

where a is the annihilation cross section, v is the relative velocity, and the bracket means 
thermal average. Here we assume that all the stop shadrons have decayed into the lightest 
one which is a neutral mesino, = tu, and similarly T°* = t*u for the anti-stops. With this 



assumption, the number density of T° is equal to that of the stop nj ll|]. We also assume 
that the number density of stop is equal to that of anti-stop. Here a and v can be simply 
expressed by cr ~ i?^ad ^'^d v ~ ^/'iT jmj . 

The timescale of the annihilation is compared with the Hubble expansion rate, 

Tann = 3i?, (7) 

where H = y/p/ {^/SMp) with Mp is the reduced Planck mass, and p is the total energy 
density, p = n'^ /SOg^^T^. We adopt effective degrees of freedom g^: = 17.25 (or 10.75), 
assuming that the QCD phase transition occurs at T = Tqcd ~ 150 MeV (or ~ 100 MeV). 
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We get the final abundance of tlie stop after the second annihilation to be 

m,Y,= 0.87 X 10--GeV ( ^V' (^Y'" ( J^Y'" (^^f\ (8) 

where the yield variable of the stop is defined by 1^ = n-jr/ s with s is the entropy density 
(=|p/T). From the temperature dependence in Eq.Q, it can be seen that 1^ is lower for a 
higher temperature after the QCD phase transition, which means that the number density 
of the stop must immediately be frozen-out just after their hadronization and annihilation at 
T = Tqcd- We see that the stop abundance can naturally be cast into the attractive range of 
10"^^ GeV ^ Bhfn^Yj ^ 10"^'^ GeV ^ for a several hundred GeV stop mass with a reasonable 
hadronic branching ratio Bh ~ C'(l)- Note that the above final abundance is lower than the 
stop abundance after the standard freeze-out, which occurred at around T ~ mj/30 11 1. 
Moreover, the decay of the thermally-produced stop would have only negligible contribution 
to the dark matter relic density I^dm^^ ^0.1 [15.] . 



IV. BIG-BANG NUCLEOSYNTHESIS AND STOP NLSP 



We consider the 2-body decay process of stop into a gravitino and a top quark. The 
emitted top quark then immediately fragments into hadrons and produce lots of high-energy 
protons and neutrons. Those emitted particles may modify the abundances of light elements 
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4l|. 



such as D, T, ^He, "^He, ^Li, ''Li and ''Be after/during the BBN epoch 

The high-energy hadrons scatter off the background protons and "^He, and induce hadronic 
shower. These processes produce energetic neutron, D, T and ^He. These non-thermally 
produced neutron and T (or ^He) then scatter off the background protons and "^He. It is 
followed by n-p (T-'^He) reactions, synthesizing D (^Li). In addition, the nonthermally- 
produced neutron also induces sequential reactions to reduce the ''Be abundance through 
a set of processes ^Be(n, p)^Li(p,^He)'^He, reducing the primordial ^Be abundance which 
means reducing the primordial ^Li abundance as well through the electron capture at a later 
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m 



time. This mechanism in non-standard BBN scenarios have been studied by 
detail. Because we are interested in a relatively short lifetime of stop ^ lO^s with its 



To be more precise it is Bh{mj~ m^)Yj: that we need to look at. However, the stop mass is much larger 
than the gravitino mass in the interesting region. 
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FIG. 2: Allowed regions in the plane of the stop mass and the gravitino mass in case of -Rhad = 
O.l/mTT. Each white region is observationally allowed for (a) ^Li, (b) ^Li and (c) D, respectively. The 
shadowed region which is left with respect to the dotted line is excluded by the experimental lower 
bound m-f- ^ 250 GeV. The left region with respect to the thick line is not allowed kinematically 
when we consider only the two-bodies decay, t ^ t + G. 



high branching ratio into hadrons {Bh ~ 1), the photodissociation processes induced by the 



emitted high-energy charged particles and photons are not important in this study [40|, 141 1. 

We would have to check whether the abundances of the other light elements also agree 
with the observational data. In particular, it would be a crucial problem if copious deuteri- 
ums are produced in this scenario. For the observational deuterium abundance, we adopt 
(D/H)obs = (2.82 ± 0.26) x 10~^ |42|. Because the stop abundance before its decay is very 
small as was discussed in the previous section, the stop decay does not significantly affect 
the "^He mass fraction and the ^He to D ratio. 

In Fig. [2] we plot allowed regions in the plane of the stop mass and the gravitino mass, 
assuming /o- = 0.1. Each white region is observationally allowed for (a) ''Li, (b) ^Li and (c) 
D, respectively. The left region with respect to the thick line is not kinematically allowed 
for the 2-body decay mode, t ^ t + G, and has stop lifetime longer than ~ 10^ s which 
we do not consider here. Note that even if we include this region our results would not be 
changed. Furthermore, there is experimental constraint which excludes mj ^ 250 GeV. 
We can understand the exclusion/ allowed regions in Fig. [2] by looking at the yield BhmjY^ 
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FIG. 3: Combined allowed region in the plane of the stop mass and the gravitino mass by using the 
constrains in Fig[2](a), (b), and (c). The name of each element is written in the respective exclusion 
regions. There is an allowed region at around Tn^ = 400 — 600 GeV and m3/2 = 2 — 10 GeV. 



and the lifetime r-j-. The analysis results in terms of these two variables is shown in Fig. 1 of 
Ref. 22]. From Eq.(8) we found that the yield of the stop grows as m^^^ as the stop mass 
increases. For small stop mass the effect on ^Li is too small, while for large mass the effect is 
too large. It is just right in between. In addition the effect also depends on the stop lifetime 
which can be seen from Fig. [1] While for ^Li, as shown in Fig. 1 of 22|, it prefers small 
yield with relatively large lifetime, or a certain range of lifetime almost independent of the 
yield. As for D, it is allowed if the yield is small or the lifetime is small. 

By combining all three constraints shown in Fig. [2] (a),(b),(c), we get the favoured region 
at around mj= 400 — 600 GeV and ^3/2 = 2 — 10 GeV which is clearly shown in Fig. [31 
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V. CONCLUSIONS 



We have studied the decays of stop NLSP to gravitino LSP and the effects on BBN. We 
found some range of stop and gravitino masses where all light element abundances agree 
with the observational data, including ^Li and ^Li which are problematic in the standard 
BBN. It is, therefore, very interesting to test this hypothesis further. The gravitino itself 
would be practically undetectable, aside from its gravitational effect, due to its very weak 
interaction. Direct detection of a dark matter particle can therefore exclude our model ^. 
We can only hope that future collider experiments, such as the upcoming LHC, would be 
able to provide some evidence on the physics beyond the standard model, and advance our 
understanding of the Early Universe. 
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